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ABSTRACT 

We consider a magnetic impurity coupled to the hybridizing and screening channels of 
a conduction band. The model is solved in the framework of poor man's scaling and 
Cardy's generalized theories. We point out that it is important to include a two-body 
hybridization if the scaling theory is to be valid for the band width larger than U. We 
map out the boundary of the Fermi-non-Fermi liquid phase transition as a function of the 
model parameters. 
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It has been known for some time that at energies much less than the Kondo temper- 
ature the ground states of simple models of magnetic impurities [1,2] are Fermi liquids 
[3,4]. Generalized versions of such models show a transition between Fermi liquid and 
non- Fermi liquid ground states [5,6,7,8]. It has recently been argued that analysis of this 
transition may lead to a better understanding of the apparently non-Fermi liquid-like 
behavior observed in the superconducting cuprates. 

Non-Fermi liquid states appear in the phase diagrams of impurity models with finite 
range interactions. A spinless impurity model has been studied by multiplicative renor- 
malization group [9] in [6] and using poor man's scaling theory [10] in [7]. For a magnetic 
impurity model, non-Fermi liquid phases have also been found [5] using Wilson's numer- 
ical renormalization group [3]. Something similar also occurs in a generalized Hubbard 
model in infinite dimensions [8]. In this paper, we adapt poor man's scaling theory and 
its extension [11] to find the phase diagram for a magnetic impurity model at finite U 
including the effect of screening channels in the particle-hole symmetric case. Until now 
scaling arguments have only been developed for the case of infinite U . We point out that 
for U less than the bandwidth, a fully renormalizable model must include a 'two-body 
hybridization process' if the renormalization group equations are to remain consistent. 
This is the process in which two electrons with opposite spin hop from the conduction 
band directly onto the local orbital. 

For finite U and in the particle-hole symmetric case, we find that there are four regions 
in the phase diagram. These are characterized by three quantities, the usual hybridization 
ti, a two-body hybridization t 2 and the spin exchange V y . The first is the non- Fermi liquid 
fixed point at which all of the quantities are irrelevant. The second region corresponds 
to what we call the L V y relevant region', where the fixed point Hamiltonian belongs to 
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the same universality class of Kondo's strong fixed point model, and is therefore, the 
Fermi liquid fixed point. The third one is the l t 2 relevant region'. We have calculated the 
exponent of the leading term of the fixed point interaction Hamiltonian and found that 
this fixed point is again the Fermi-liquid fixed point. We find that the region where all 
°f 9\i 9y an d #2 are relevant can be divided into three sections. Two of them have the 
same fixed point as that of the ^-relevant and t 2 - r elevant regions. In the other section 
the system behaves as 'free orbital Fermi liquid' [2]. 

Our model Hamiltonian is 

H = Yl e kc{ a iC ka i + e d n d + Un d ^n dl 

k>0,a,l 

+ ti J2(cl d a + h.c.) + t 2 J2{c< a0 c<_ a0 d a d- a + h.c.) (1) 

cr a 

+ Vic ] al c a i$ a ,d a > + ^ Vxic^Caidlda + ^ V y ^ al C- a iS_ a d a . 

a,a' ,1 <r,l <r,l 

The first four terms in (1) are the usual Anderson model while Vi term takes into account 
the finite range interactions between the local l d orbital' and conduction electrons (/ = 
0,...,Nf). The V X) i and V Vi i terms describe the spin exchange interaction and are taken 
to be zero for / = 1, ...,Nf for convenience. The t 2 term [12] describes processes in which 
two I = 0-channel conduction electrons with opposite spin hop onto the local orbital. To 
renormalize the model in the framework of poor man's scaling theory it is important that 
this term as well as spin exchange terms are included. 

Following the terminology of [5,6], we call the 0-channel the hybridizing channel and 
the others the screening channels. We deal with the screening channels by bosonization. 
The spin of the screening channels is not important and we treat the screening channel as 
spinless for simplicity. Let bki be the bosonic operators corresponding to for / = 1, Nf 
[13]. The effective Hamiltonian then is H = H + where 
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H o = e kcl a c ka + e d n d + U n d ^n di 

k>0,a 

k 

+ v o cU*<i>d*> + Vx Yl c l c * d l d ° + - b li b kh (2) 

a,a> cr fc,/>0,<7 P 

Hj = h E( At 4^ + h.c.) + t 2 J2(^ 2 ct c -* d ad-a + h.c.) + V y J2 clc_ a dl a d a , 

where the index denoting the / = channel has been suppressed. The electron dispersion 
is taken to be = (k — kp)/p with p = (/i-yp)" 1 being the density of states at the Fermi 
surface. The operator A is given by A = exp{-E k pVi/VkL{b kl -bl)}. The Hilbert space 
of H can be projected onto four subspaces characterized by the four possible impurity 
electron states, \a >— |0 >, \a > and |3 >= | |!>- Each term in the interaction Hi plays 
the role of a dipole operator causing transition among those subspaces. 

The Hamiltonian in (2) can be treated using Haldane's familiar procedure [14] by 
writing the partition function in terms of a sum over histories of the impurity. Expanding 
the partition function in Hj and labelling the (imaginary) times so that < T\ < ... < 
r n < (3, the partition function can eventually be written as a sum over all possible histories 
of the local degrees of freedom which fluctuate between the 4 local states. This gives 

Z = Y1 J2 I — ... / — exp(-5 , [ri, ...,r n ; a u ...,a n }), 

n a u ...,a n J T J T 

(3) 

S[r,a]=Y E CiC + >V"^ h ^ )+ ^^^' 

i<j a=cr,l i i 

where a±, a n (a n+ i = a±) and t±, r n label a Feynman trajectory, r is the ultraviolet 
cut-off. With a simple shift of the ground state energy, E a can be chosen so that J2 a E a = 
0. We take E = E 3 = — E a = —e d r/2 for the particale-hole symmetric case (e^ = —U /2). 
g a p are coupling constants, with g 0a = g a3 = g± = hr, g a - a = g y = V y r and g 03 = g 2 = 
t 2 r. 
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Equation (3) may be thought of as describing a one-dimensional four-component 
plasma of kinks carrying 'charges' q a a p and 'fugacities' g a/3 [11,8]. E a is a 'magnetic field'. 
The 'charges' are given by q% a = ((1 - S x /n)Sx a - So/n^i/n), q a Ga , = ((1 - 5 x /ir)(5 Xa/ - 
5 Xa ),0), q a a3 = ((1 - 5 x /7r)5 x ^ - <S /MiA) and g a 3 = ((1 - 5 x /n) - 25 /n, 25,/tt). 
g a/ g = — ^q,. The phase shifts are 5 — 2 tan" 1 7rpVo/2, 5 X = 2 tan -1 npV x /2 and 5; = 
2 tan" 1 npVi/2 [17]. The 'charges' obey the relations (? Q/ g + <jig 7 = g a7 , which means 
that the model can be regarded as a special case of the general one-dimensional model 
with 1/r 2 interaction considered by Cardy [11] and can be renormalized by poor man's 
scaling theory. The Coulomb term in (3) can be rewritten as a one-dimensional spin 
chain model with interaction K(aj, aj)r 2 / (r, — r,) 2 , where K(a,/3) = — 1/2 X)<x(<22/?) 2 > 
ie. X(0, cr) = X(cr, 3) = -70, K(a,a') = -j x (l - 5 a ,o') and K(0,3) = ~f x - 2 7 o with 
K(a, (3) = K((3, a) and K(a, a) = 0. Here 7o = (1 - 5 x /n - 5 /n) 2 + (5 /tt) 2 + £,( V*") 2 
and 7 X = (1 — 5 x /ir) 2 . It is worth noting that if we consider the model with vanishing t 2 , 
the relations q a p + qp-y = q a -y are not satisfied so that the model can not be mapped to 
the Cardy's model. 

Cardy has pointed out that the critical behavior of these kinds of model can be 
discussed using a generalized poor man's scaling theory. Si and Kotliar have extended 
Cardy's theory to the case in which small 'magnetic fields' E a are included [15] as Hal- 
dane has also done [14]. For the particle-hole symmetric case of the present model, the 
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renormalization group equations are given by 

"' rt = (i - \io)gi - gig y e- E ° - gig2e~ Eo , 



d\nr v 2 

</f/2 = (l-2 7o + lx )g 2 - 2gfe- L - 



amr , , N 

(4) 



= -4 lx (gle- E °+g 2 y e- E °), 
= e d r - 2gj(e- E ° - e~ E °) + 2gle~ E « - 2g^ 



d\n.T 
de d r 
dlnr 



We make some comments on these equations. The first three equations of (4) are exact 
for any value of e d r while the last three are restricted to the case where the impurity level 
e d is much less than the band width 1/r [14]. ( The assumption of a large band width is 
within the spirit of poor man's scaling theory.) The last equation reflects the variation of 
e d with 1/r: 

" -2^-*° - e-*0 + 2g\e- E « - 2 9 y E ° . (5) 



d lnr 

The requirement E a = which we enforce throughout the scaling process, ensures 
that the separation of the field e d r renormalization and the free energy renormalization 
is unambiguously determined. 

We have assumed that g 0a = g a3 , K(0,a) = K(—a,3) and K(0,3) = AK(0,a) — 
K(a, —a) in the present model. The scaling must preserve these equalities if the model is 
to remain consistent. This only happens in the particle- hole symmetric case. Away from 
the particle-hole symmetry, the scaling generates different flows for g 0(T and g a3 , K(0, a) 
and K(—a, 3), etc. 

As we have already said, if the model can be mapped to the special case of Cardy's 



model, t 2 can not vanish. Furthermore, from (4), we can see that even if we start with 
vanishing g 2 and g y , their absolute values all increase at a rate proportional to g ± . This 
means that the flows calculated assuming that g y and g 2 vanish exactly are not correct 
renormalization flows. To describe a fully renormalizable model when the band width 
is larger than U, the Hamiltonian must include these two-body hybridization and spin 
exchange terms. 

Since our renormalization group is perturbative in its treatment of g a p, the renormal- 
ization of 70 and r y x can be neglected in the first three equations of (4). This allows us to 
draw out the phase diagram in 7 — ^ x space (See Figure 1), which can be divided into 
four regions: 

(i) For 7q > 2, 7 X > 1 and 27 — j x > 1, all g a p are irrelevant and renormalize to zero. 
There exist weak coupling fixed points g* af3 = 0. The fixed point Hamiltonian is similar to 
the multi-channel X-ray edge problem [16]. The system exhibits a power-law decay of the 
correlation function with a non-universal exponent. This is the no n- Fermi liquid phase. 
Since 70 = (1 — 5 /ir — 5 x /ir) 2 + (<5 /vr) 2 + I];(5//7r) 2 , it is clear that it is the existence 
of the screening channels which allows 70 to exceed 2 and leads to the non-Fermi liquid 
phase [7]. 

(ii) For 70 > 2, 270 — 7^ > 1 but 7 X < 1, <7i and g 2 are irrelevant and renormalize to 
zero. g y is relevant. We focus on the fixed point Hamiltonian. Choosing the parameters 
V y = Jj_, V = 2J Z and V x = —J z with J z < (this ensures 7^ < 1), we see that the fixed 

point Hamiltonian H* is just the single-channel Kondo Hamiltonian: 

k 

H* = e k c L c ka + e d n d + Un d] n dl + ~ b li b ki 

k>0,a k>0,l=l P (g) 

+ J -±(Sts- d +h.c.) + J z S* c s* d . 
Here S c and s d are the spin operators of hybridizing and localized electrons respectively. 
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This g y ( JjJ-relevant region is controlled by the strong-coupling fixed point of the Kondo 
problem. Furthermore, the values of the interaction parameters V and V x are not impor- 
tant as long as they are within the parameter regime of the g y relevant region. Hence, the 
fixed point corresponding to this g y relevant phase is the Fermi liquid fixed point. Near 

the fixed point equation (5) reduces to, r — — - = —2gf , which implies that U increases. 

a lnr y 

This is consistent with the results for symmetric Anderson model: when the impurity level 
width r ~ t\ is much less than U, the model is equivalent to the Kondo model [2,18]. We 
therefore call the ^-region the Kondo region. 

(iii) For 70 > 2, 7^ > 3 but 270 — 7 X < 1, <7i and g y are irrelevant and g 2 is relevant. 
This parameter regime corresponds to the empty and doubly-occupied states being favored 
over singly-occupied states. We call this the G 2 -region. The fixed point Hamiltonian 
H* = Hq + Hj where can be read off from (2) by replacing all parameters by their 
fixed point values. Hj has its leading term R = (gl/r) J2a( c a c -adad- a + h.c.) . The matrix 
elements of R in the fixed point basis can be evaluated by analogy with the X-ray edge 
problem [19]. As was done in [5], we denote by \n d = > the eigenstates of |n d =o an d 
by \n d = 2 > the eigenstates of H*\ nd=2 . Then 

<n d = 0\R\n d = 2 >~ r a . (7) 

The anomalous exponent a — (1 — 270 + 7x)/2. If a < the operator R is irrelevant, 
and the fixed point corresponds to the first kind of fixed point we have discussed. In 
the present parameter regime, a > and the operator R is relevant. The fixed point is 
regarded as a Fermi liquid one [5]. The line a = is marginal. (We note in passing that 
if V x = 0, and hence j x — 1, our result appears to recover that of [5]. However, as the 
parameter region in which (7) holds requires 7^ > 3, this may be just a coincidence.) 

(iv) For 7 < 2, all coupling constants are relevant. We divide this parameter region 
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into three sections (Figure 1). Although all coupling constants are relevant, (4) shows that 
9y ^ 92 m 92 ^ 9y m II ari( i 9y ~ 92 m HI- According to (5), the Hubbard interaction 
U increases and so singly-occupied states are favored in I which is in the Kondo strongly 
coupling phase. In II, U decreases and the empty and doubly-occupied states are favored. 
This suggests that in 77 the system is the same FL phase as the G 2 -region. In ///, U 
remains close to its initial small value and |0 >, \a > and |3 > are mixed. The impurity 
level width T(~ tf) ^> U in III implies that the system is in the 'free orbital Fermi liquid 
phase' [2]. 

We may summarise what we have learnt about the phase diagram as follows. There 
are two kinds of fixed point in the phase diagram: Fermi liquid and non-Fermi liquid. 
The part of the phase diagram controlled by Fermi liquid fixed point can be divided into 
three regions according to the behaviour expected at finite temperature. The three types 
of behaviour are Kondo strong-coupling, free-obital and what we have called 'GV- 

It is interesting to compare our results with those of the numerical renormalization 
group reported by Krishna-murthy, Wilson and Wilkins [2]. The case they studied cor- 
responds to a point in our phase diagram: 7 = 7 X = 1. For this case, we reproduce 
qualitatively their results from the renormalization group equations (4). They found that 
the system goes to a Fermi liquid fixed point with two regions: which they call free or- 
bital region (this corresponds to our region /// in Fig. 1) and local moment region (I 
in Fig. 1). Because they assume that t 2 is always zero, it is clear that they would not 
find the region corresponding to our region //. The 'GV' and 'NFL-' regions in Fig. 1 
are physically relevant. If there is a direct Coulomb interaction between the conduction 
electrons and the impurity electron (i.e. Vq ^ and V x ^ 0), the system could be driven 
into the //-region even if the two-body hybridization vanishes initially. 
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In conclusion, we have derived that the scaling laws for a magnetic impurity model 
including the 'hybridization' of an up and a down spin electron hopping onto the local 
orbital. The renormalization group analysis shows that there is a Fermi-non-Fermi liquid 
transition when the strength of the local interaction is varied. For the case of finite U 
which we have considered, the model has to be particle-hole symmetric in order to preserve 
the consistency of the renormalization group equations derived from the generalized poor 
man's scaling theory. 

The asymmetric model can be discussed only in the infinite U limit, when the doubly- 
occupied states are completely suppressed leaving only local states |0 > and \a >. This 
is then just a special case of the spin-iV + 1 model with added screening channels and 
has been discussed by Si and Kotliar [15]. Results of a numerical renormalization group 
calculation were also reported in [5]. At finite U, the normal assumption that the usual 
hybridization remains particle-hole symmetric after renormalization is valid only for the 
particle-hole symmetric case. Results of more general treatment will published separately 
[20]. 
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Fig.l The phase diagram in 7 — 7 X space. The thicklines are phase boundaries. The 
thin lines divide the Fermi liquid phase into the different regions characterised by 
different behaviour at finite temperature.. 



13 









G 2 


4 








lx 




II 




2 




/ 111 


NFL 




111 / 


I 


KONDO 



2 4 



70 
FIGURE 1 



14 



